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1. Introduction 

Molecules are associated with the 'cold' Universe. Historically astronomers first studied 
the spectra of hot, bright stars full of atomic emission and absorption lines. Except for the 
nearby Sun, cooler stars with molecular lines are fainter and more difficult to observe. As 
the surface temperature drops, the star reddens and the peak of the emission curve shifts 
to the infrared. (The Wien displacement law for a blackbody, AT = 2898 /um K, gives a 
maximum at 1 /jm for 2898 K. [1]) The total radiated power also varies approximately as 
the fourth power of temperature (from the Stefan-Boltzmann law in so far as stellar 
emission is approximated by a blackbody. The spectroscopic study of 'cool' sources with 
surface temperatures in the range of about 500 — 4000 K is therefore best carried out 
with large modern 10 m class telescopes to gather as many photons as possible in the 
near infrared and infrared spectral domains. 

Molecules are in fact ubiquitous in the Universe and are found in diffuse clouds, dark 
clouds, giant molecular clouds, planetary nebulae, circumstellar envelopes, stellar at- 
mospheres, brown dwarf atmospheres, planetary atmospheres, comets, supernovae and 
galaxies. Molecules in dark clouds, circumstellar envelopes, and star-forming regions are 
usually studied by recording pure rotational spectra using radio telescopes because these 
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long wavelengths are scattered less than infrared and optical wavelengths in these dust- 
obscured sources. Temperatures are very low, typically less than 50 K. There are already 
a number of good reviews on interstellar clouds and related objects from an astrochemical 
perspective 

aaaa 

and they will not be covered here. Instead this review will focus on 
visible and infrared observations of molecules in cool stars, brown dwarfs and extrasolar 
planets using a mixture of spectroscopic and astronomical citations. The distinction be- 
tween molecules found in stellar atmospheres and those in molecular clouds is not always 
clear because many stars are enshrouded by clouds of molecules and dust, which can 
then appear in the stellar spectra. For example, Mira variable stars have prominent hot 
water overtone absorption bands that originate from at least two pulsing circumstellar 
shells of water vapour 0]. 

Modelling of the chemical composition of stellar atmospheres is also different from that 
of dark clouds. The very low temperatures of dark clouds mean that only reactions with 
very small energy barriers are feasible (e.g. ion-molecule reactions) and the source is 
not in chemical equilibrium 0]. The composition of a dark cloud is therefore calculated 
by assuming an initial composition and solving the coupled differential rate equations 
as is typical in chemical kinetics. In contrast, the higher temperatures and pressures 
in stellar atmospheres mean that the system has largely attained chemical equilibrium. 
The composition of a shell in a stellar atmosphere can therefore be predicted given an 
elemental composition, temperature (T), pressure and Gibbs free energy of formation 
(AGf(T)) of all possible compounds as a function of temperature [81]. In this case, the 
coupled algebraic equations that describe the equilibria need to be solved simultaneously 
to determine the composition, as is customary in chemical thermodynamics. 



2. Astronomy Background 

Big Bang nucleosynthesis created only the elements H, He and a trace of Li some 13.7 
billion years ago [9|, [lOj; all other elements are formed in stars 0,0. After the Big 
Bang, the early Universe expanded rapidly and cooled. After a few hundred million years, 
gravitational attraction produced the first stars, the 'Population III' class of objects 



13 ]. Little is known about this first generation of stars because they have not been 



detected yet with certainty; a few molecules such as LiH are predicted to form in the 



early Universe LJ] . Models predict that population III stars were massive (more than one 
hundred solar masses) and their intense UV emission ionised the surrounding interstellar 
medium. Massive stars burn the H and He fuel rapidly to form the elements up to Fe in 
the periodic table by nuclear fusion reactions. The Fe nucleus is the most stable in terms 
of binding energy per nucleon so fusion reactions are unable to make any heavier elements. 
After a few million years of life Population III stars began to collapse as the nuclear fuel 
was exhausted, and they exploded as supernovae. Nuclear reactions during supernova 
explosions release a rapid burst of neutrons (r-process) that are absorbed by nuclei and 
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make elements heavier than iron. Some heavy elements are also produced by nuclear 
reactions in stellar cores from the slow neutrons (s-process) released in fusion reactions. 
Supernova explosions and stellar winds therefore return matter to the interstellar medium 
enriched in 'metals'. In astronomy all elements other than H and He are termed metals. 

The Population II class of stars were created by the gravitational collapse of matter in 
the interstellar clouds that formed from the metal-enriched ejecta of Population III stars. 
Population II stars are still depleted in metals compared to younger stars like our Sun 
that are part of Population I. The elemental abundance of a star is a crucial property and 
our own Sun is taken as a standard reference. Astronomers have defined the metallicity 
of a star as the abundance ratio of iron to hydrogen as compared to the same ratio in 
the Sun, i.e. the metallicity [Fe/H] is defined as 



[Fe/H] = log 10 — = log 10 -rj^ - log 10 C 1 ) 



(N Fe /N u ) Sun J ^ V J star V / Su „ 

so a metallicity of corresponds to the same relative Fe abundance as the Sun. Fe is cho- 
sen for measurement convenience, and metallicities of —1 and —2 correspond to relative 
Fe depletions by factors of 10 and 100, respectively, compared to the Sun. Population 
II stars have metallicities more negative than —1 and Population III stars should have 
essentially no metals and metallicites less than —6, i.e. Fe depleted by more than a factor 
of a million, are expected. A few low mass stars with metallicites as small as —5.3 have 



been detected, but their evolutionary histories are not known with certainty [15]. 



Normal stars are classified primarily by colour 16]. This so called 'main sequence' 
contains the usual stellar types (O B A F G K M) classified by temperature from about 
60 000 K (blue, early O-type stars) to 2400 K (red, late M-type stars). Each letter is 
divided into 10 subdivisions (0 — 9) so at 5800 K our Sun is a G2 dwarf. 'Early' M's 
(or other stellar types) refer to low numbers and 'late' M's to cooler higher numbers. 
The classification of stars is two dimensional, based both on surface temperature (e.g. 
G2 for 5800 K) and 'size' (luminosity). In practice there are three main categories for 
luminosity: supergiant (I), giant (III) and dwarf (V), although five Roman numerals (I to 
V) are used. Luminosity (with 1 L = 3.845 x 10 26 W equal to the luminosity of the Sun) 
is used for convenience instead of the more difficult to determine mass (in units of solar 
masses, 1 M & = 1.989 x 10 30 kg). Dwarfs are the largest category and our Sun is thus 
classified as a G2 V star. Stellar interiors have temperatures of millions of degrees (e.g. 15 
million K for our Sun), but they are not visible from the exterior because of the 'opacity' 
of stellar atmospheres. Implicit in the stellar main sequence is the correlation between 
mass and surface temperature: massive stars burn their nuclear fuel more rapidly and 
have higher surface temperatures. 

In practice stellar classification is based on the presence or absence of key atomic and 
molecular features in the spectrum [It]], which is sometimes called the spectral energy 
distribution. A-type stars have strong hydrogen lines and the coolest stars have strong 
near infrared bands of electronic transitions of diatomic TiO and VO. Recently, objects 
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cooler and smaller than M-type stars have been discovered and called brown dwarfs 18]. 
The classification system was therefore extended with two new spectral types, L and T; 
some early L's are stars, and late L's and T's are brown dwarfs. L-type objects have 
temperatures of about 2400 — 1300 K and are characterised by the presence of near 
infrared electronic transitions of the metal hydrides FeH and CrH, and the fading of VO 
and TiO bands. T-type brown dwarfs have temperatures in the range of about 1300 — 700 
K and have strong absorption from overtone vibration-rotation bands of methane and 
water. Finally a new 'Y' category of very cool brown dwarfs has been predicted with 
temperatures below about 600-700 K and characterised by strong ammonia overtone 
absorption bands [l9T | . 

Brown dwarfs are not stars because they are not massive enough to burn hydrogen in 
their cores, although nuclear fusion of other nuclei such as deuterium and lithium is still 
possible. Effective surface temperatures are still somewhat uncertain, but the new field 
of observation and modelling of sub-stellar objects is progressing rapidly. 

The term 'dwarf is used in several different contexts in stellar astronomy. For example, 
the term subdwarf has little to do with luminosity but refers to a low abundance of heavy 
elements 20]. A subdwarf has a metallicity of at least —1 and an extreme subdwarf has 



a metallicity of —2. In fact the lack of metals in cool subdwarfs means that the opacity 
of the stellar atmosphere is decreased so that subdwarfs can be more luminous in the 
visible than dwarfs (i.e. subdwarfs are not necessarily 'sub-luminous'). The terms dwarf, 
brown dwarf and subdwarf therefore refer to different kinds of objects. 

Not all stars lie on the 'main sequence' discussed above. Clearly very young stars — 
young stellar objects or YSOs — that have just been born in dark interstellar clouds do 
not lie on the main sequence. YSOs will evolve until they are main sequence stars and 
then spend most of their life on the main sequence [jj. Similarly when low-mass stars 
similar to our Sun convert most of their H nuclear fuel into He, they evolve off the main 



sequence and become red giants 2l|, |22j. Red giants burn He in their cores to make C, 



and their outer shells expand and cool. The red giant phase is followed by a complex 
AGB (asymptotic giant branch) phase where the core consists of mainly C and O with 
energy coming from the fusion of He and/or H in thin surrounding shells 21]. During the 
AGB phase strong stellar winds return much of the matter to the interstellar medium 
and the object evolves towards a planetary nebula. AGB stars often possess circumstellar 
shells of dust and molecules which envelop the star and render it nearly invisible in the 
optical region. In a planetary nebula the circumstellar shell has become detached from 
the hot central stellar core. A planetary nebula is a spectacular short-lived phase (a 
few thousand years) in which the surrounding cloud of matter glows in the visible from 
the intense UV excitation of the very hot central star. Ultimately the cloud of matter 
dissipates and the central star becomes a white dwarf. 
High-mass stars with stellar masses greater than about 8 M & will explode as Type II 

r ii 

supernovae, unless mass- loss is sufficiently rapid [23]. In essence, high-mass stars consume 
their nuclear fuel until a sufficiently large (> 1.3 M@, the Chandrasekhar limit) inert 
Fe/Ni core remains in which nuclear fusion is no longer possible. The heat production 
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of the core is insufficient to support the outer layers and gravitational collapse ensues, 
rapidly releasing an enormous quantity of energy. The outer layers of the star disintegrate 
and ultimately the core ends up as a neutron star or a black hole. 

The ratio of carbon to oxygen abundance, C/O, is a crucial parameter. Normal stars 
on the main sequence have C/O less than 1 and when molecules begin to form CO is 
always present because of its strong triple bond. The formation of CO ties up most of the 
carbon and an oxygen-rich chemistry results. If C/O is greater than 1 then the formation 
of CO ties up all of the oxygen and a very different carbon-rich chemistry is possible. Red 
giant and AGB stars produce carbon from helium so the composition of stars that are no 
longer on the main sequence becomes increasingly carbon rich. When the abundance of 
carbon and oxygen are equal then S-type stars (see below) result [24|] and when carbon 
is more abundant than oxygen then carbon stars (C-type stars) form 25f] . Carbon stars 
are therefore red giant or AGB stars with C/O > 1. 

Giant planets such as Jupiter (1 Mj = 0.955 x 10~ 3 M & = 318 M©), which has a 
surface temperature less than 200 K, are believed to fall on the same main sequence as 
brown dwarfs. Hundreds of extrasolar planets have been discovered primarily by the very 
small periodic Doppler shifts they induce in the absorption lines of the spectrum of the 



parent star 26|, |27J] • A few extrasolar planets have been detected by the small decrease 



in luminosity observed when they transit in front of the star. There is a selection effect 
that preferentially favours the discovery of 'hot Jupiters' because massive planets with 
short periods are easier to find. Extrasolar planets are therefore closely related to brown 
dwarfs and will be included in the set of cool objects covered in this review. 



3. The Sun and Sunspots 



As our nearest and most important star, the Sun holds a special position in astronomy 
and its study is part of a separate sub-discipline. The Sun has an effective surface tem- 
perature of 5800 K, which dissociates all polyatomic molecules. Lines of neutral atoms 
and atomic ions dominate the solar spectrum as first observed by Fraunhofer in the 
visible region. Fraunhofer had no explanation for the dark absorption lines in the solar 
spectrum and labelled them with letters starting with A in the near infrared at 760 nm. 
Apart from the A- and B-lines (which are the forbidden A- and B-bands of terrestrial 
O2, i.e. the 0-0 and 1-0 vibrational bands of the b x Sg — X 3 £~ electronic transition), all 
of the Fraunhofer lines are due to atoms; for example, the C-line at 656 nm is the first 
H Balmer line (H a ) and the D-line contains the two Na D-lines (589.0 and 589.6 nm). 

A handful of diatomic molecules contribute weakly to the near infrared, visible and 
near UV solar spectra 0,0: MgH A 2 II-X 2 £ + 0,0, C 2 Swan bands (d 3 II g -a 3 II u ) 



3J, CH B 2 IT - X 2 n and A 2 A - X 2 n [3J, CN B 2 £+ - X 2 S+ (Violet System) [3J and 
A 2 n-X 2 S+ (Red System) [351]. Solar spectra in this region have been recorded with the 
1-m Fourier transform spectrometer (FTS) at the McMath-Pierce Solar Telescope on 
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Kitt Peak and published in a series of solar atlases by L. Wallace and co-workers 26l. |291] ; 
there is also an earlier high resolution solar flux atlas 361 ] . 

In the infrared, the solar spectra are more interesting as they contain prominent molec- 
ular bands. Excellent ground-based spectra are again available from the Kitt Peak FTS 

r~\ I — T I — 1 

covering the 1850-9000 cm -1 37J, |38[] and 460- 630 cm -1 regions [39(]. In addition, ob- 



servations from Jungfraujoch in Switzerland (a higher, drier site than Kitt Peak) extend 
to longer wavelengths [40J] and cover 250-630 cm -1 . At short wavelengths the CN Red 
System and the C2 Phillips System (A 1 !^ — X 1 ^^) |4lf] can be seen, but the spectra are 
dominated by the strong first overtone and fundamental vibration-rotation bands of CO 
near 2.5 /xm and 5 ^m 42J. The Meinel bands (OH vibration-rotation bands) near 3 urn 
are also strong and the OH pure rotational lines are prominent at long wavelengths 43 ] . 
The vibration-rotation bands of CH can be seen at 3 (im 33j, [44]. 

In the infrared, excellent high resolution solar spectra have also been recorded from 
low earth orbit by the ATMOS and ACE FTSs. The ATMOS instrument produced 
two solar atlases covering 625-4800 cm" 1 at 0.01 cm resolution [45|, |46|] and the 



ACE atlas covers 750-4400 cm -1 at 0.02 cm^ 1 resolution 47J, |48] with a higher signal- 
noise ratio and no telluric interference from residual gases in the instrument, Figure 1 
( http: / / www.ace.uwaterloo.ca/solaratlas.html [ ) . The satellite spectra also clearly show 
vibration-rotation (49| and pure rotational lines of NH. The combination of laboratory 
and solar spectra have led to a considerable improvement in the spectroscopic constants 
of OH 43], NH 49], and CH 33] in their electronic ground states. 

Sunspots are cooler than 5800 K — hence their apparent contrasting black colour. Large 
black sunspots (Figure 2) can be as cool as 3300 K, equivalent to an M4 star, but 
with a vertical magnetic field of more than (X3 T. There are a series of Kitt Peak at- 
lases for sunspots that cover 470-1233 cm -1 
5l|, 4000-8640 cm" 1 Q, 8900-15 050 cm- 1 



39(] ^970-8640 cm" 



50], 1925-3480 cm" 1 
53j, and 15 000-23 000 cm" 1 [54]. Pa- 



per copies are available from L. Wallace (wallace@noao.edu) and electronic copies from 



f tp : / / nsokp . nso . edu/pub/ atlas / . 



All of the molecules seen in the photosphere intensify in the spectra of sunspots, with 



the exception of C2, which fades and is not obvious. The weak MgH B' 2 S + — X 2 £ 



r 2 V + 



31 



551 ] transition makes its appearance in the near infrared, along with the A x n — X 1 !] 



transition of A1H 
CaH in the red 
abundance ratios 



56J] near 425 nm and the B 2 £+ -X 2 £+ and A 2 n-X 2 S+ transitions of 
57 1 . The MgH lines can be used to determine 24 Mg: 25 Mg: 26 Mg isotopic 
551 ] . The lines of several electronic transitions of TiO [60] become 



extremely dense in the visible with the 7-bands (A 3 $— X 3 A) 58] and 5-bands (b 1 n— a 1 A) 



591 ] in the red and near infrared. The Sun and sunspots are remarkable sources for 



spectroscopy at 3000-6000 K, and the two TiO papers by Ram et al. [581 . [59J nicely 
demonstrate how laboratory and solar data (Figure 3) can be combined to give improved 
spectroscopic constants. 

As in M-type stellar spectra the FeH F 4 A — X 4 A transition 6l|, [62] near 1 //m and 
the E 4 n — A 4 n transition 63J near 1.58 yum are present in sunspots, but with a large 
Zeeman effect. The FeH F-X 0-0 band near 990 nm is called the Wing-Ford band in the 
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Figure 1. The ACE solar spectrum with the intensity scale normalised to 1. Strong CO absorption 
is seen near 2000 and 4000 cm -1 and OH pure rotational lines near 1000 cm -1 are clear. 



astronomical literature (see below). Interestingly the sunspot spectra were combined with 
laboratory observations to derive rotational ^-values for levels involved in some of the 
low- J lines in the Wing-Ford band (641]. These ^-values are useful in measuring magnetic 
fields in cool objects. For hotter objects such as the Sun and sunspots the magnetic field 
strength is obtained from magnetically-sensitive atomic lines (i.e. lines between levels 
with large (/-values) either by the observation of Zeeman splittings for large fields or 
from the polarization of the lines for weaker fields. In cooler stars and brown dwarfs 
there are no useful atomic lines, so it has been proposed to use the Zeeman effect on the 
lines in the Wing-Ford band instead. 
The infrared spectra of sunspots display a great profusion of lines: in addition to the 



CO, OH, NH and CH molecules seen in the photosphere, SiO [65J, HF [66J], HC1 [67]], and 
H2O 68|,[69( appear. The fundamental and first overtone bands of SiO cover a substantial 



part of the 10 (im and 5 //m regions 65J], and H2O has an amazing density of up to 50 
lines/cm -1 starting with the pure rotational region at the long wavelength limit and 
extending to nearly 7000 cm -1 [70j] in the near infrared. 

An important advance in the spectroscopy of hot water occurred as a result of the 
detection of water vapour in sunspots. At 5800 K, the Sun's photosphere is too hot for 
water to exist, but by 3900 K the concentration of OH and H2O are equal. A large 
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Figure 2. An image of a sunspot as seen with the Vacuum Tower Telescope of the National Solar 
Observatory at Kitt Peak. 



number of unassigned lines were noticed in two Kitt Peak sunspot atlases 39, 50]. It was 
suspected that these lines were due to hot water but the available laboratory data were 
inadequate to confirm this. New laboratory infrared emission spectra of hot water at 
1800 K were therefore recorded. Comparison of the laboratory emission spectra of H2O 
and the sunspot absorption spectra identified most of the unassigned sunspot lines as 
H2O lines, but only a small fraction of the new water lines could be assigned quantum 
numbers. This work proved that there was 'water on the Sun' but the line density was 
so high that the spectrum seemed 'unassignable' by conventional techniques 68] . 

We began a collaboration with the theoreticians Polyansky and Tennyson to apply 
more sophisticated approaches to the problem [g^]. Through variational calculations of 
the energy levels using a high quality ab initio potential energy surface, Polyansky et 
al. 69J were able to assign most of the strong lines. The assignment method relies on 



the smooth variation of errors and confirmation of tentative assignments by combination 
differences. 

More recently, we have recorded a laboratory spectrum of water vapour at 3000 K by 
using an oxy-acetylene torch as a source. The water emission in the 500-13 000 cnf 1 
spectral region was recorded with a high-resolution FTS 7l[ . Work on this 3000 K spec- 
trum of H2O has just finished with the publication of the last paper in the series covering 
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Figure 3. The 0-0 <5-band of TiO as (a) seen in the laboratory in emission as (b) seen in a sunspot 
in absorption and (c) absent in the solar photosphere [59| . Reproduced by permission of the A AS. 



the near infrared part [72|]. This spectrum matches closely the conditions encountered 
in the photospheres of late M dwarfs and sunspots. A paper on 'monodromy', that is, 
the change in the bending energy level pattern observed when the highly-excited water 
molecule starts to sample linear geometries, has also been published [73|]. In this paper, 
we identify vibration-rotation energy levels for the highly excited bending states 8f2 and 
9^2 above the barrier to linearity for the first time. In our analysis we included sunspot 
spectral data. 

An important application of hot water work is the calculation of molecular water 
opacities for various cool objects. Reliable potential and dipole surfaces are needed to 
compute the millions of transitions required to reproduce low-resolution infrared spectra 
of M-type stars and brown dwarfs. The latest water line list of Barber et al. 7J| contains 
more than 500 million lines and is recommended for the simulation of the spectral energy 
distribution of cool stars and brown dwarfs. 

Recently the solar abundances of C, N, O, and Ne have been revised downwards by 
nearly a factor of two based on the intensity of atomic and molecular lines through the 
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use of sophisticated 3-dimensional radiative hydrodynamic model atmospheres 75|, |76| . 



771 . 1 781 ] . The metallicity of the Sun was decreased from the previously accepted value of 



about 0.0245 down to 0.0165 [78J] (for the Sun, the metallicity is defined as the mass 
fraction of the elements that are not H and He) . 

The granulation (outside of the sunspot in Figure 2) seen in solar images is caused by 
the upward motion of matter in the hotter centre of a gas cell ('granule') and downward 
motion at the cooler edges. The complex fluid motion in the solar atmosphere leads to 
unusual asymmetric line shapes and to line shifts. The line width is larger than expected 
based on the Doppler formula ll|] because of mass motion, traditionally attributed to 
'micro' and 'macro' turbulence [79J]. In fact modern radiation-hydrodynamic models are 
capable of reproducing the observed line shapes [80J, but it is not clear yet whether 
the new abundances are correct 
models of the Sun 



811 ] . The lower abundances seem to be inconsistent with 



771 ] used in the relatively new field of helioseismology 82J that observes 



periodic Doppler shifts of solar lines. These Doppler shifts are due to the propagation of 
sound waves in the Sun and have a typical period of 5 minutes. 



4. K and M Stars 



The Sun and particularly sunspots cover much of the temperature range spanned by 
K (4000-5200 K) and M (2400-3700 K) stars. Perhaps the best published spectrum is 



that of Arcturus (a-Bootis) , one of the brightest stars in the northern sky [83j, [84j, |85 ] . 
Arcturus is a K1.5 Illpe red giant star with a surface temperature of 4320 K (pe stands 
for 'peculiar emission' because of the presence of emission lines in addition to the usual 
photospheric absorption features). The infrared and near infrared atlas (0.9 — 5.3 /im) of 
Arcturus was recorded with the FTS that was associated with the 4-m Mayall Telescope 
831 ] at Kitt Peak, the visible atlas (373-930 nm) with the 0.9-m coude feed telescope 
and the coude spectrometer at Kitt Peak [84| and the UV atlas (115-380 nm) recorded 
mainly with the STIS instrument on the Hubble Space Telescope [85J. As expected, based 
on the solar atlases, the same molecular features are present with the exce ptio n of C2, 



and the Lyman and Werner bands of H2 86] appear in emission in the UV 

M stars are particularly important because they are the most numerous type of star 
in our Milky Way galaxy, with 70% by number and 40% by mass 87]]. Because they are 
so faint, it is only relatively recently that a full range of M dwarf spectra have become 
available 0, 89]. The visible and near infrared spectra of M-type stars are dominated 
by TiO absorption bands (Figure 4) as discussed above for sunspots. In Figure 4, the 
terms 'active' and 'inactive' refer to magnetic activity which induces atomic emission 
lines such as Balmer H a at 656 nm, and is associated with stellar rotation [ixj. Notice 
that astronomers have retained the old atomic notation with Na I referring to the neutral 
Na atom and Ca II to C& + . 
The TiO 0-0 A 3 $ - X 3 A band head appears near 705.4 nm in the Ml spectra (3700 
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K) and has nearly become saturated by M6 (2750 K) with weaker TiO features such as 
the 0-1 band at 758.9 nm now present [9J]. By M6 the B 4 II - X 4 £~ 1-0 and 0-0 bands 



92|, |93J of VO are present near 733 nm and 785 nm (Figure 4) and strengthen for late 
M objects (Figure 5) Q. As in sunspots MgH (B' 2 S+ - X 2 £+), CaH (B 2 £+ - X 2 £+ 
and A 2 n - X 2 £+ at 635 and 690 nm, Figures 4 and 5 ) and FeH (F 4 A - X 4 A) are all 
prominent. Perhaps the most interesting new feature (Figure 4) in M spectra is ascribed 



to the CaOH A 2 n - X 2 S+ transition at 625 nm [94l . 1951 ] . In the infrared the overtone 
spectra of water appear more strongly than in sunspots for the late M's (Figure 6). 

In fact it was in late M dwarfs that Wing and Ford 96| were the first to detect a mys- 
terious band near 991 nm at low resolution. The Wing-Ford band was later detected in 
S-type stars and in sunspots at higher spectral resolution [97|. Nordh et al. [98J| identified 
the Wing-Ford band as the 0-0 band of a FeH electronic transition by comparison with 
an unassigned laboratory spectrum that had a band head at 989.6 nm. The 1-0, 2-0, 2-1 
and 0-1 bands of this transition can all be identified in sunspot spectra. The spectra of 
the Wing-Ford band and the other bands are very irregular in appearance because of 
many overlapping branches with perturbed lines. Quantum chemistry predicts a large 
number of low-lying electronic states and their mutual interaction causes perturbations 
that appear as irregular patterns of lines. These 'many-line spectra' with no obvious pat- 
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Figure 5. Red and near infrared M9, L3 and L8 dwarf spectra 102]. Reproduced by permission 
of the A AS. 



terns at first sight are characteristic of the electronic spectra of transition metal hydrides. 
After heroic efforts, the near infrared transition of FeH was shown to be a F 4 A-X 4 A 
transition and seven bands with v', v" < 2 were assigned by Phillips et al. 6l| These 
bands were recorded in emission using the FTS at the National Solar Observatory at 
Kitt Peak using a carbon tube furnace source at about 2300° C. More recently the line 
intensities for the F 4 A-X 4 A transition were obtained through ab initio calculation [g^]. 

The bands of CaH are particularly useful in identifying M subdwarfs that have reduced 
metal abundances compared to our Sun 99]. TiO and VO which characterise M-stars 
contain two heavy elements ('metals' in the astronomical sense), while CaH has only a 
single heavy element. If the metal abundances are low, then the CaH bands strengthen 
relative to TiO, and the ratio of CaH to TiO band intensities can be used to identify 
subdwarfs [991 ]. 
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Figure 6. Near infrared M and L dwarf spectra [1031 ] . Reproduced by permission of the AAS 



5. Brown Dwarfs 



The discovery of the first brown dwarf Gl 229B by Oppenheimer et al. lOOf ] in 
1995 was a momentous event and now hundreds of substellar objects are known 
(http://DwarfArchives.org/). Brown dwarfs are actually defined by mass because ob- 
jects with less mass than about 0.075 M = 79 Mj (1 M = 1047 Mj, mass of Jupiter) 
are unable to fuse hydrogen in their cores but can still burn deuterium and lithium 



101] ]. There is a similar deuterium burning limit at about 13 Mj and the lithium limit 



is about 63 Mj. Cool objects in the L and T classes, however, are classified primarily 
by effective surface temperature, which controls the appearance of the spectra, and not 
by their difficult-to-determine mass. While there is certainly a correlation between mass 
and surface temperature, the relationship is not simple because of the different possible 
ages of the objects. Low mass objects burn hydrogen, deuterium, lithium, etc. slowly and 
can therefore have lifetimes of billions of years, and indeed some lifetimes are calculated 
to be comparable to the age of the Universe. It is therefore possible for a small young 
brown dwarf to have the same surface temperature (and spectral class) as a larger old 
low-mass star that has been slowly cooling. 
The L class of objects is defined by the weakening of the VO and TiO metal oxide bands 
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Figure 7. Near infrared T dwarf spectra [l07|. Reproduced by permission of the AAS. 



and the strengthening of the FeH and CrH metal hydride bands (Figure 5) [lsl . Eq3 - The 
CrH A 6 S + -X 6 S + transition with a 0-0 band he ad at 861 nm first appears in late M's 
(Figure 5) and is probably not present in sunspots 104J . 1 105f | . For later L's, the absorption 
bands of CrH and FeH fade and are replaced by infrared absorption of H2O and CH4 
overtone and fundamental vibration-rotation bands (Figures 6 and 7) 18|, Il06l ] . 

The determination of surface temperatures of brown dwarfs is more complicated than 
for stars. In part the problem is due to their faintness and to difficulties observing in 
the infrared with ground-based telescopes. As Figure 7 illustrates, the spectral energy 
distributions no longer resemble blackbody curves — molecular absorption has chopped 
the spectral energy distributions into pieces. Surface temperatures therefore need to 
be determined using simulations and these models are still relatively primitive with 
incomplete molecular opacities. The transition from M to L objects is particularly difficult 
to model because of the condensation of gases to form dust; the L to T transition is easier 
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to handle because the dust has settled out and has disappeared. In L-type objects particle 
scattering needs to be included in the simulation of spectral energy distributions in order 
to match the observed extinction. Existing effective surface temperatures are ~2500 K 
for L0 and ~1500 K for L9 while TO start at ~1500 K and descends to ~800 K by T8 



18fi. 



Brown dwarf masses are also difficult to measure. Evolutionary models predict that 
while L's could be stars or brown dwarfs, T's are always brown dwarfs. One possible 
independent way to determine mass is from the surface gravity; astronomers are often 
able to determine the temperature, metallicity (composition) and gravity from high res- 
olution stellar spectra. Gravity affects the observed spectra by changing the thickness 
and structure of the photosphere, which changes the appearance, f or examp le, of the 
wings of lines. The first high resolution spectra are starting to appear 10g, [109J, but one 
problem is that all known brown dwarfs are rotating rapidly with velocities in excess of 
20 km/s, which corresponds to a spectral line broadening due to the Doppler effect of 
about 1 cm -1 near 1 fim 110]. 

The importance of high spectral resolution is illustrated for TiH. Opacity calculations 
predict that TiH shou ld be particularly a bund ant for late M-type subdwarfs |lll| . In 2003 
Burgasser et al. [112l ] and Lepine et al. [1131 ] classified two objects (2MASS 0532+8246 
and LSR 1610-0040, respectively) as metal-deficient L-type subdwarfs. In the Burgasser 
et al. paper 1121 ] there is an 'unassigned hydride' band at 960 nm and another unassigned 
band at 940 nm. Both of these features are overlapped by a water overtone band in the 
Earth's atmosphere and hot water absorption in the object s the mselves. In a later paper 
on a third L-type subdwarf (2MASS 1626+3925) Burgasser [1141] suggests that the fe ature 
at 940 nm might be due to the A 4 $-X 4 $ transition of TiH 1 151 ] . Reiners and Basri 116 ] 
then recorded higher resolution spectra (resolving power, A/AA, of 31000) of two of the 
three L-subdwarfs. They conclude that LSR 1610-0040 is not an L-type subdwarf but 
a peculiar slightly metal-deficient M-type dwarf with strong CaH bands and relatively 
weak CrH and FeH features. The feature near 960 nm is due to a group of Ti lines and 
there is no sign of the TiH lines near 940 nm in the high resolution spectrum of the L 
subdwarf 2MASS 0532+8246. 

For T-dwarfs the main contributors to the opacity are water, methane and ammonia. 
Cushing et al. 1171 ] recorded a nice sequence of moderate resolution spectra of M, L 
and T dwarfs in the 0.6 to 4.1 /im region with ground-based telescopes. These spectra 
show beautiful methane features at 1.65 /jm (21/3, first overtone of the C-H stretch, 
Figu re 8) that begin to saturate at 3.4 ^m (1/3, C-H stretch fundamental) for a T5 dwarf 



llgj. Near 2.3 //m the spectra also show the CO first overtone bands for L's fade as the 



temperature decreases towards the T's and the dominant carbon-containing compound 
becomes methane (Figure 9). Water is present in all the L's and strengthens towards the 
T's (Figure 9). More recentl y spe ctra to longer wavelengths have become available from 
the Spitzer Space Telescope 1191 ] and they clearly show the umbrella mode of ammonia 
1^2 at 933 cm -1 (10.7 //m) in T dwarfs (Figure 10). 



Unlike the situation for water 



74 ]. reliable molecular opacities for hot methane and 
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Figure 8. Spectra of L dwarfs showing the FeH E-A transition and of T dwarfs with the methane 
C-H first overtone band [l 1 Tj ] . Reproduced by permission of the A AS. 



hot ammonia are not available. Hot emission spectra of methane at 800, 1000 and 1273 K 
were recorded to match the range of temperatures encountered in T dwarfs 1 181 ] . These 
spectra are not assigned so it is di fficult to use them except at the temperatures at which 
they were recorded. Cushing et al. 1171 ] have used our hot laboratory spectra to identify 
features due to methane in T dwarf spectra. 

The best existing database for methane and ammonia is HITRAN 2008 1201 ] . which 
is aimed at the calculation of infrared transmission of the Earth's atmosphere at 300 
K. HITRAN is missing the hot bands needed for high temperature sources and is also 
rather incomplete in the near infrared. New laboratory spectra for hot methane and 
ammonia are clearly needed and some work is underway in our research group. There is 
also steady progress on the calculation of the methane spectrum using ab initio methods 



1211 ] . but the current predictions are far from spectroscopic accuracy; similar theoretical 
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Figure 9. Near infrared L and T dwarf spectra showing CO first overtone bands 1171 ]. Reproduced 
by permission of the AAS. 



work has been started on ammonia 1221 ] . Hot ammonia is particularly significant for the 
next co oler spectral class — Y brown dwarfs — that have surface temperatures less than 
the T's 



123]. 



6. AGB, S and C Stars 



The spectra of S- and C-type stars are complicated because of the variety of molecular 
species that can form as the C/O abundance ratio increases as stars evolve off the main 
sequence and start to burn He. S-type stars have C/O of about 1 while C-type stars have 
C/O greater than 1 and are defined empirically by the presence of certain molecular 
bands in their spectra. These stars have surface temperatures that are roughly parallel 
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to those of M stars and are believed follow the evolutionary sequ ence M-MS-S-SC-C, 
allowing for the possibility of intermediate MS an d SC stages 1241 ] . Values for C/O are 



about 0.5 for M, 0.6 for MS, 0.8 for S and 1 for SC [125l |. As discussed in the introduction, 
at the end of a star's life on the main sequence stellar evolution proceeds through the 
red giant, AGB, planetary nebula and white dwarf stages. The empirical M-MS-S-SC-C 
sequence is therefore a consequence of stellar evolution and corresponds to the red giant 
and AGB part of stellar demise. 

S-type stars are traditionally classi fied by the pres ence of ZrO and YO bands in addition 
to the usual TiO bands of M-stars 



124 11251 . 



126]. The heavy elements Y and Zr are 



beyond the Fe abundance peak and are formed by s-process neutron reactions that are 
characteristic of stars that have left the main sequence. The classification notation used 
is Sx,y or Sx/y, in which x is the usual to 9 for surface temperature and y is a number 
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from 1 to 7 based on t he re lative intensity of TiO and ZrO bands (which is a proxy for the 
increasing C/O ratio) 126 ]. A number of visible t rans i tions of Z rO (/3 system, c 3 II — a 3 A; 



7 system, b 3 $ — a 3 A and B 1 !! — X 1 E + t ransition 
and A 2 I1 — X 2 S + transitions 



129 



124 



122, W) and YO (B 2 S+-X 2 E+ 



130l |) are seen. With much of the oxygen tied up in 



CO, a number of unusual su lfur-contai ning mole cules are able to form in S-type stars 
including TiS 



131 



13l ll33j and ZrS [134 



135|l observed in the near infrared in, for 



130])- 



example, R Andromedae (R And, Figure 11 

A high resol ution FTS spectrum of the S-type star R Andromedae was recorded by 
Ridgway et al. 1371 ] in the 2400-2800 cm -1 region. Vibration-r otat i on li nes of OH, NH, 



CH and HC1 fundamental bands and the first overtone of CS 
and later the fundamental band of SH 



140 



13a . Il39| ] were identified 



1411 ] was noted. To shorter wavelengths CO 



and the Red System of CN are very strong in v Cygni, an S-type Mira variable that has a 
regular pulsing Stella, atmosphere Q. Low resolution spectra from the Infrared Spaee 
Observatory (ISO) 1251 ] suggest that HCN can be seen near 3 microns 1431 . [144]], and 
recen tly SiO and SiS were identified in Spitzer Space Telescope spectra near 10 microns 



145]. 
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There are two types of S stars: intrinsic and extrinsic that can be sepa rated by the 
presence of the unstable element technetium (Tc) in the photosphere 146]. Radioactive 
"Tc (the isotope formed by neutron reactions in AGB stars) has a half-life of 2.1 x 10 5 
years and does not occur naturally on Earth. Tc, however, is an s-process element that 
can be identified through its atomic lines in intrinsic S stars (i.e. genuine AGB stars). 
Extrinsic S stars do not have Tc because it has all decayed away. It turns out that all 
extrinsic S stars are binary with a white dwarf as a companion. Apparently the heavy 
s-process elements such as Y and Zr in extrinsic S stars were transferred through mass 
loss during the companion's AGB phase. 

Carbon stars have perhaps the richest of all stellar spectra and this has caused consid- 
erable confusion in their classification. Traditionally C stars are identified by the strong 
visible abs orption o f the C2 Swan system, CN Red and Violet Systems and the CH A-X 



transition [147L Il48l . 1 1491 ] . Carbon stars can also display the Merr il l — S a nfor d bands due 
to the A 1 B2 — X 1 Ai transition of SiC2 (origin band near 497.7 nm 1501 ]) and the comet 
system of C3 (A 1 ^ — X 1 !]^ transition near 405.2 nm 151 ]). The simple approach of 
classifying carbon stars by their surface temperature and the amount of carbon (as mea- 
sured by the strength of the C2 Swan bands) is not very successful because of the wide 
range of objects that have C/O greater than 1. The surface temperatures of C stars are 
notoriously difficult to determine. The current system divides C stars into three main 
sub-classes R, N and CH (denoted as C-R, C-N and C-H), and then further categorises 
each of the se by temperature, strength of C2, peculiar carbon isotope ratios and lumi- 
nosity 1481 ] . The R type is hotter and has no Ba lines from the s-process, while the 
cooler N-type has enhanced heavy elements with strong blue absorption. The CH type 
has strong CH absorption and the J-subtype also has strongly enhanced 13 C abundances, 
with 12 C/ 13 C ratios as high as 2 as compared to the solar value of 92. Isotopic abun- 
dances show great deal of variation in astronomy and often provide important clues to 
the evolutionary history of the object. CH carbon stars are believed to be part of binary 
systems with white dwarfs (like extrinsic S st ars) and their anomalous abundances are 
due to mass transfer from the other star 



1521 ] . 



The near infrared spectra of carbon stars are d omi nated [1531 ] by the CN Red System, 
the Phillips and Ballik-Ramsay (b 3 £~ — a 3 n u 1541 ]) systems of C2 and CO overtone 
bands. The vibration-rotation bands of HCN and C2H2 appear in the 3 fim region and 
to slightly longer wavelengths CH, NH and CO fundamentals and the CS firs t overton e 
can be seen in both ISO spectra [1551 ] and at higher resolution in FTS spectra [1371 . Il56l ]. 



In the ISO spectra [1551 ] the fundamental band of CS and the first overtone of SiS are 
detected near 7 (im plus t he v i bending mode of HCN at 14 /jm and the ^5 bending 
mode of C2H2 at 13 .5 /im. [1571 ] . The det ection of C3 has also been suggested using low 



1581 ] and Spitzer spectra 



159, 



Ir- 



resolution ISO 

It is also possible to have oxygen-rich AGB stars and their infrared spectra typically 
show OH, H2O, SiO, CO, CO2, SO2, HC1 and HF, as expe c ted f rom their composition 
and the surrounding cloud of molecules from mass loss 



161 



162.1163]. 
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7. Extrasolar Planets 

The first extrasolar planet — a 'hot Jupi ter 1 orbiting around the star 51 Pegasi — was 



discovere d by Mayer and Queloz in 1995 [1641 ]. the same year the first brown dwarf was 



observed 100]. The parent star is classified as G2 V and is very similar to the Sun, but 
the 51 Pegasi b planet is very strange indeed. (Planets are named using the parent star 
and lower case letters b, c, d, and so forth in order of discovery.) 51 Peg b has a mass 
of about 0.5 Mj, but is located at 0.05 AU (1 AU is the mean Earth-Sun distance of 
1.5 x 10 11 m) and has an orbital period of 4.2 days. In essence it would be as if Jupiter 
were to be transported so it was 7 times closer to the Sun than Mercury in our Solar 
System. The composition of 51 Peg b is expected to be similar to that of Jupiter but 
with a surface temperature of about 1300 K rather than about 125 K for Jupiter. 

Almost all of the hundreds of extrasolar planets (or exoplanets, see 
|http:/ / exoplanet.eu/| have been detected by observing periodic Doppler shifts in 
the spectrum of the parent star. A star and a planet orbit about a common centre of 
mass and this motion causes the parent star to periodically move away and toward the 
observer. For example, the Sun- Jupiter system would have a periodic Doppler shift of 
13 m/s for lines in the Sun [27|. Detection sensitivities are currently better than 3 m/s, 
which corresponds to a frequency shift of about 5 MHz for the Na D line. Interestingly, 
most of the Doppler shift observations view the star through a temperature-stabilised 
iodine cell (as a frequency reference) and record a high resolution visible spectrum. It 
is not possible to determine the inclination angle of the orbital plane from the Doppler 
effect, but these observations give the orbital radius, eccentricity of the orbit and a lower 
limit for the mass of the planet (as long as the stellar mass is known). The Doppler 
shifts are largest for large planets close to the parent star so there is a selection effect 
that favours the detection of hot Jupiters. Given a few technical improvements the 
Doppler shift method should be able to detect Earth-like planets. 

For extrasolar planets, it is possible that the orbital inclination (the angle between the 
Earth-star vector and the perpendicular to the star-planet orbital plane) is close to 90° . 
The planet can then transit the star and be detected by the small decrease in light as the 
planet passes in front of the star. The hot Jupiter HD 209458b was the first transiting 



planet observed with a dip in light intensity of about 1.5% [1651 ] . This type of precision 
photometry is best done from orbit and the Corot (http://smsc.cnes.fr/COROT/) and 
Kepler missions (http://kepler.nasa.gov/) are currently searching for Earth-like planets 
by their transits. The combination of Doppler and transit data allows for a rather com- 
plete description of the planet including orbit inclination, planetary mass, radius and 
density. 

The technique of 'transit spectroscopy' allows the detection of molecular absorptions 
in exoplanets. The idea is simple: the effective size of a planet with an atmosphere 
depends on the wavelength because the planet appears bigger at a strongly absorbing 
wavelength than at a wavelength for which the atmosphere is transparent. By measuring 
very precisely the size of the dips in the stellar radiation during transit as a function 
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of wavelength, a spectrum of the planet is obtained. Swain et al. 1661 ] observed the 
transiting exoplanet HD 189733b with the Hubble Space Telescope using the NICMOS 
camera in the 1.4-2.5 (im range. HD 189733b is a hot-Jupiter planet in a 2.2 day orbit 
around a bright K2 V primary star. A model atmosphere of mainly molecular hydrogen 
with a water relative abundance of 5 x 10 -4 and a methane abundance of 5 x 10 -5 was 
needed to match the observed wavelength dependence of t he tr ansit dips. They used the 
BT2 [74J linelist for water and a combination of HITRAN 1201 ] and Nassar and Bernath 



[1181 ] for methane. 

More recently the 'emergent flux' from the same planet was directly dete cted 1657 ] 
with a similar instrument configuration as used for transit spectroscopy 166]. The NIC- 
MOS instrument on Hubble was used again for imaging spectroscopy with a grism (a 
grating ruled onto a prism inserted in the field of view of an infrared camera) that gave 
a resolving power of 40 in the 1.4-2.5 /jm range. In essence Swain et al. manipulate the 
spectra recorded just before the planet disappears behind the star (i.e. star plus planet) 
and spectra recorded during occulation with the planet behind the star (i.e. star only) to 
obtain the emission spectrum of the planet alone. The planetary flux contained modula- 
tions that required the addition of H2O, CO and CO2 to their model atmosphere. The 
interpretation of such spectra are not straight forward because the molecular bands could 
appear in absorption or emission or both, depending on the details of the atmospheric 
pressure-temperature height profile and the composition of the different atmospheric lay- 
ers. The non-detection of methane is explained as being due to the different atmospheric 
regions sampled by emission spectroscopy (dayside emission from higher pressure and 
temperature regions) as compared to transit spectroscopy (sampling of lower pressure 
cooler outer regions in the terminator between night and day). 



8. Molecular Opacities 

The spectra of cool objects are rich in molecular absorption features, which presents a 
major difficulty in the calculation of model atmospheres needed to simulate the spectral 
energy distribution functions. Compared to atoms, molecules have many more energy lev- 
els and hundreds of millions of absorption lines are needed for molecular opacities. It is 
not feasible to obtain millions of line parameters from laboratory observations alone: they 
must be provided by theory. Yet the purely theoretical calculation of molecular opacities 
is also not satisfactory because, except for molecular hydrogen, ab initio quantum chem- 
istry is not yet sufficiently accurate. In other words, laboratory measurements alone do 
not provide sufficient lines while purely theoretical calculations do provide enough lines, 
but not with enough accuracy. The solution to this problem is to extend laboratory 
measurements with theory to provide the required molecular opacities. In essence, the 
strongest lines seen at high resolution come from laboratory data, but the millions of 
weaker lines that provide a kind of 'continuum' absorption come largely from theory. 
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The calculation of molecular opacities has been recently reviewed by Sharp and Burrows 



168]. 



Our approach to the calculation of molecular opacities for the electronic transitions of 
metal hydrides such as FeH illustrates the contributions chemical physics can make in 
astronomy. The simulation of these spectral energy distributions require spectroscopic 
line parameters (positions, intensities and pressure-broadening parameters) to be used 
in a model atmosphere 1691 ] . The generation of the required molecular opacities is not 
an easy task for molecules such as CrH, FeH and TiH. The assignment of perturbed 
laboratory spectra are difficult and require considerable patience using combination dif- 
ferences. The use of just laboratory data, however, is not satisfactory for astronomical 
purposes because even for metal hydrides tens of thousands of lines contribute to the 
observed stellar and substellar spectra. The vibrational and rotational levels generally 
need to be extended to higher v's and J's than observed in the laboratory. This can 
be done empirically or by using the predictions of ab initio quantum chemistry. Purely 
ab initio opacities, however, are not satisfactory because the calculated line positions 
are not accurate enough. A judicious combination of experiment and theory works best 
to cover a wide range of v and J. Minor isotopes and satellite branches should also be 
included, usually from theory, because when absorption from the main lines saturate, 
the weaker satellite and minor isotopologue lines become important. 

Line intensities are also required and the best approach is to measure the radiative 
lifetimes for a few lines by, for example, laser spectroscopy. The transition dipole moment 
function can be computed by ab initio quantum chemistry, but typical errors for state-of- 
the-art calculations are 10-25%. For e xamp le, our prediction of the radiative lifetime of 
the A 6 E + v=0 level of CrH was 0.75 /jls [1041 ] and a subsequent experimental measurement 



1701 ] yielded a value of 0.939 /is, some 25% larger. Perhaps the most reliable approach 



is to compute the ab initio transition dipole moment function and then scale it with 
an experimental measurement to match the observations. Theory is thus being used 
to extend an experimental value to cover a larger range of vibrational levels. Finally, 
pressure-broadening parameters can also be determined by experiment or theory, but 
they are diffic ult to measure or compute so generally some reasonable values are just 
assumed 



1051 ] . 



The first step in obtaining molecular o pacit ies for CrH, FeH and TiH is to use the 
observed lines for the A 6 £+-X 6 £+ (CrH F 4 A-X 4 A (FeH [41]), A 4 $-X 4 $ (TiH 



1151 ]) and B 4 r-X 4 <3? (TiH 1711 ]) transitions in fits to obtain the usual spectroscopic 



constants for each vibrational level. These constants are obtained by deweighting the 
perturbed lines, which is particularly difficult for the heavily perturbed FeH bands. The 
next step was to carry out ab initio calculations of molecular properties using multiref- 
ere nce config uration interaction (MRCI) method with the MOLPRO suite of programs 



62, 



104, 



llll ]. Large atomic basis sets were used and scalar relativistic effects were in- 



cluded. In the initial multireference part of the calculation, the valence electrons were 
distributed among the valence orbitals ('active space') to generate a large number of 
Slater determinants that were treated equally in a CASSCF (Complete Active Space 
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Self Consistent Field) calculation. Additional electronic correlation was then included by 
configuration interaction (CI). These extensive state-of-the-art calculations are required 
to obtain reliable molecule properties for transition metal hydrides. 

The calculated potential energy points for each electronic state were used to calcu- 
late vibrational wavefunctions, vibrational energies and vibrationally-averaged rotational 
constants (B v ) by solution of t he on e dimensional radial Schrddinger equation. This was 



particularly important for TiH [11 11 ] because only the 0-0 vibrational bands are available 
for the A-X and B-X electronic transitions. For TiH, therefore, the vibrational constants 
were taken from ab initio calculations for each electronic state. The vibration-rotation 
energy levels were cal culat ed (including the effects of spin-orbit and spin-spin coupling) 



for FeH [62| and TiH 111| for v's from t o 4 or 5 and rotational levels typically up to 
N of about 50. In our first effort on CrH 1051 ] . we covered only up to v=3 and about 
N=40, which should be extended to higher v and N. These tables of energy levels are 
adjusted to contain the experimental term values derived directly from the observations 
if they are available, otherwise calculated values are used. In this way, the experimental 
line positions are reproduced by the energy levels. 

In the next step, all possible transitions are computed from the tables of energy levels, 
even the weak satellite lines and weak bands not seen in the experiments. Each calculated 
line requires an intensity so an Einstein A{j value was calculated. The starting point was 
the transition dipole moment for each electronic transition calculated from the electronic 
wavefunctions. The vibrational wavefunctions were then used to calculate the Einstein 
A value for each vibrational band. The Einstein A value (in s _1 ) for each line was 
then calculated using the formula A = A V ' V »HLF/(2J' + 1), in which HLF is a Honl- 
London factor (for rotation) and A v i v n is the Einstein A value for the v' — v" band. The 
Honl-London factors were obtained from the Hund's case (a) expressions (which had to 
be derived for the A- A, 4 $- 4 $ and 4 r- $ cases) and the rotational wavefunctions 
obtained during the calculation of the energy levels. The line strength [l| of each line 
is most convenient in the form of an integrated cross-section J adv in units of cm 2 s™ 1 
molecule -1 . The integration over frequency is just to eliminate the line shape function. 
The integrated cross-section for a line is calculated from the Einstein A by 

f adv = — \—A(2J' + 1) exp(-E"/kT)(l - exp(WfcT)) (2) 
J 8iri> z Q 

with D in cm -1 ; Q is the internal partition function. As the partition function appears in 
the line strength function, we also re-evaluated the thermochemistry for each molecule 
using the best estimates of the dissociation energies. We computed the partition function 
using estimated spectroscopic constants for all of the low-lying spin components and 
electronic states. The correct evaluation of the partition function is a frequent problem 
in astronomical applications and the simple rigid rotor, harmonic oscillator expressions 
that are approximately valid at 300 K, are rarely reliable at high temperature. 
The new molecular opacities and thermochemistry for CrH, FeH and TiH have been 
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incorporated into a spectral synthesis code by Burrows 62|, Il05l . A comparison of 
the observed spectrum of the L5 brown dwarf 2MASS-1507 and the spectral simulation 
(including FeH and CrH) was carried out 62J. The agreement was reasonable, considering 
that L dwarfs are hard to model because of the presence of dust in their at mospheres . 



More recent comparisons with high resolution spectra are not so favourable [I08l 11091 ] 
and it is clear that many of the line positions and intensities will need to be adjusted 
based on experimental observations. 



9. Conclusions 

Astronomy of the cool Universe is the study of molecules. As illustrated in this review, 
the tools of modern chemical physics can be used to great advantage in astrophysics. 
What is also evident is that astronomy in turn provides some unique environments for 
molecules and for the study of molecular processes that are rarely encountered in the 
laboratory. A remarkable early example was the observation of only the predissociated 
lines of the A 1 II — X 1 !]" 1 " t rans ition of A1H in emission in the spectrum of the MS- 



type variable star x Cygni [1721 ] . These A1H lines are not seen in laboratory emission 
spectrum because of tunnelling through a potential barrier in the A state, but appear 
through inverse predissociation in the stellar atmosphere as Al and H atoms recombine. 
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